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Abstract 

Lattice power-counting is extended to QCD with staggered fermions. As 
preparation, the difficulties encountered by Reisz's original formulation of the 
lattice power-counting theorem are illustrated. One of the assumptions that 
is used in his proof does not hold for staggered fermions, as was pointed out 
long ago by Liischer. Finally, I generalize the power-counting theorem, and the 
methods of Reisz's proof, such that the difficulties posed by staggered fermions 
are overcome. 
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1 Motivation and summary 



Lattice QCD with improved staggered fermions (SFs), or, Kogut-Susskind fermions 
[1-3], has recently enjoyed pubhcity for its abihty to correctly reproduce many as- 
pects of hadronic physics with reasonable accuracy [4,5]. However, SFs have some 
notable properties. For instance, SFs do not entirely overcome the fermion doubling 
problem. Rather, they reduce the number of continuum modes from 16 to 4. (A 
further reduction to 2 modes is possible, by projecting quarks and antiquarks to odd 
and even sublattices resp. [3].) These 4 modes are referred to as tastes, to distinguish 
them from the Nf flavors in the continuum theory. To estimate the fermion measure 
of Nf continuum flavors, one takes the power Nf/A of the fermion determinant in 
the definition of the functional integral. I will not address the attendant controversy, 
but rather another technical question: lattice power-counting for staggered fermions. 
Here again, fermion doubling creates difficulties, as will be discussed at some length 
in this article, and as was pointed out some years ago by Liischer [6] . 

To better understand perturbative renormalization of SFs it is of course useful 
to have a lattice power-counting theorem. By way of analogy, renormalizability of 
SU{N) Yang- Mills coupled to Wilson fermions has been proven some time ago by 
Reisz [7]. This result was based on his earher work on BPHZ-like renormalization 
theory on the lattice [8,9]. That work rested crucially on his lattice power-counting 
theorem [10,11]. (This literature is rather mathematical; more accessible reviews are 
those by Reisz [12,13] and Liischer [6].) Reisz 's lattice power-counting theorem was a 
significant achievement because on the lattice Feynman integrands are trigonometric 
rather than rational functions of momenta; this can lead to results that differ from 
those of the continuum in important ways.^ 

It is often stated that no power-counting theorem exists for SFs; for example in 
Refs. [15-18]. However, it is also widely believed that the theory of SFs coupled 
to Yang-Mills (denoted here SF-QCD) yields the right quantum continuum limit in 
perturbation theory. That is to say, the lattice perturbation series can be renormalized 
and matched to a continuum renormalization scheme at every order in the gauge 
coupling g. This conclusion is supported by an analysis of the types of non-irrelevant 
operators that are allowed by the symmetries of SF-QCD. One finds that all such 
operators are already present at tree-level. (See for example [19] and refs. therein.) 

^An amusing example occurs, for instance, in naive discretizations of supersymmetric quantum 
mechanics [14]. 
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That is, from a Wilsonian point of view one concludes that SF-QCD is in the same 
universahty class as continuum QCD. It is reasonable to believe that by an adjustment 
of the bare parameters of the lattice action, one can arbitrarily adjust the coefficients 
of all non-irrelevant operators in the infrared, in order to obtain the desired theory. 

The belief that SF-QCD is renormalizable also follows from a consideration of 
powers of the lattice spacing a that arise in vertices and propagators of the theory, 
and how they appear in loop diagrams, an early example being [3]. In fact, for 1-loop 
diagrams, it is easy to power-count by partitioning the loop integration domain in a 
sensible way and estimating the integrand and measure for each of those domains. 
But this is nothing other than a limited version lattice power-counting. So, in fact, 
a version of power-counting already exists, though it is not as general as we would 
like. In actuality, this sort of partitioning is exactly what is done in Reisz's proof 
of his lattice power-counting theorem. However, the complexities that occur at high 
orders — where the number of domains increases factorially — are best addressed by 
a more sophisticated mathematical approach, just as in the continuum proofs of 
Weinberg [20] or Hahn and Zimmermann [21]. It is this sort of general method of 
power-counting that is aimed at in the present study. 

I now summarize the remainder of this article: 

• In ^ I briefly review two well-known formulations of the SF-QCD action, and 
the corresponding Feynman rules. 

• In ^ I review the conditions for the Reisz power-counting theorem. I also 
remind the reader of the lattice UV degree (of divergence) that is defined in 
Reisz's theorem. 

• In 21 the conditions of the Reisz theorem are examined for the two formulations 
of SF-QCD thate were described in ^ It is shown that in both cases the condi- 
tions of Reisz's theorem are violated. I explain the essential, basis-independent 
reason for this failure. 

• In ^ I generalize Reisz's theorem and methods of proof in such a way that 
lattice power-counting can be applied to staggered fermions. 

• In ^ I conclude with a summary and discussion of further issues that could be 
explored. 

Various appendices are included for details that would detract from the main discus- 
sion, but that are essential to the proof: 
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• In discuss resolutions of identity that are used in a domain decomposition 
for the loop momenta integration. Both the one used by Reisz, and a generalized 
one that is applied in ^are given. 

• In ^^I discuss a simplification of the Feynman rules that is very useful in the 
momentum-space taste basis (MSTB). 

• In [23 it is shown that the domain of internal momenta can be extended in a 
useful way in the MSTB. 

2 Bases 

2.1 The 1-component basis 

The gauge covariant SF action with link fields U^{r) is just [1-3] 

+ X ma^x(r)x(r), a^(r) ^ {-y^+-+-^-K (2.1) 

rGZ4 

Color indices are suppressed, r is a site index, and a is the lattice spacing. I refer to 
this as the 1-component basis (ICB). Under the lattice translation 0(r) — > 0(r + s), 
with (f) = U^, X, Xi the action is only invariant for even shifts s G 2Z^. Thus what I 
will call the Kdhler-Dirac lattice^ consists of elementary cells that are hypercubes of 
length 6 = 2a on each side, denoted 2aZ^. 

The free fermion propagator [U^ = 1) has 16 lattice poles; i.e., minimal eigenval- 
ues of the (Euclidean) SF Dirac operator. The 16-fold degeneracy corresponds to 4 
continuum Dirac fermions. In momentum space, the additional poles lie at edges of 
the first Brillouin zone 

Ba = {-7i/a,n/a]\ (2.2) 

More specifically, the poles lie at the sites of the lattice (7r/a)Z^ that is reciprocal to 
the Kahler-Dirac lattice 2aZ^. 

^This is the lattice generated by the basis vectors 2/i. I.e., those corresponding to translations that 
leave the lattice action invariant. It is the lattice through which free SFs are related to Kahlcr-Dirac 
fermions [22-26]. 
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Perturbation theory is, as usual, defined by expansion of 

U^{r) = expiagA^{r) (2.3) 
in powers of g. Fourier transforms are defined witli conventions: 

m = ^ <P^{X,X,A,}, (2.4) 

wliere Ba is defined in fl2.2p . I have chosen to make k dimensionful, since it agrees 
with the conventions of Reisz. It will prove useful below to periodically extend the 
fields: 

4){k + {27r/a)z) = 0(A;), V z e Z\ (2.5) 

The form of Af^{k) follows the convention of [19] but differs from the convention 
of [27,28] by a factor of e*"'^''/^. The choice that is made here gives Feynman rules 
that are manifestly periodic on Ba since (j){k + (27r/a)z>) = 0(A;) Vz/ G {1,...,4}, 
and (f) G {x, x, ^^}- In the conventions of [27,28], slightly more effort must be 
expired to demonstrate periodicity of numerators of Feynman diagrams on Ba, since 
individual vertex factors lack this property. 27r/a-periodicity of the numerator of 
Feynman integrands is an important assumption in Reisz's proof. See for example 
the Reisz-Liischer conditions VI and CI in ^ 

For instance, in the conventions of [27,28] the 0{g) gluon-quark vertex is propor- 
tional to cos{pij_a + ^k^a), with p incoming momentum on the x line and k momentum 
on the incoming gluon line. Suppose is a loop momentum and we want to check 
the 27r/a-periodicity in condition VI. Under A; — > + (27r/a)/i the vertex reverses 
sign. This corresponds to v4^(fc + (vr/a)/!) = — y4^(fc), due to the additional factor of 
g2afc^/2 ^Y^e Fourier transform of [27,28]. However, one finds that the sign is always 
cancelled in some other part of the diagram that also involves the loop momentum 
k. This of course must be true, due to the equivalence with the formulation that I 
choose, where the Feynman rules themselves enjoy periodicity on Ba- 

The propagators and leading boson-fermion vertices are given, for instance, in 
Table 1 of [27], apart from the factor e"*'^^'^/^ for each incoming y4^(fc). Let 5a denote 
the 27r/a-periodic (5-function. Also define /2 = J2u=i ^- Thus 

1 = (0,0,0,0), 2 = (1,0,0,0), 3 = (1,1, 0,0), 4 = (1,1, 1,0). (2.6) 

These are useful because in (I2.1|l we can write a^{r) = expiirfi ■ r, which makes the 
Fourier transform easy to compute. 
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Ignoring ghosts and pure YM vertices,^ the Feynman rules are: 



-I -1 



(2.7) 



(2.8) 



Qj 



(2.9) 



Tj^ are generators of the gauge group in the quark representation R; i,j are color 
indices; A is an IR regulating mass for the gluon. p, q are incoming momenta on the 
X, X lines respectively, whereas k, £ are incoming momenta on lines. Note that 
momentum conservation is only mod it /a where the fermions are concerned (i.e., one 
finds (7r/a)/i inside the 5a-functions) , due to the fact that the Kahler-Dirac lattice is 
2aZ^ (in physical units), which has for a reciprocal lattice ('7r/a)Z^. 

2.2 Momentum space taste basis 

Here I discuss the momentum space taste basis (MSTB) that was originally introduced 
in [3]. I present the results in conventions that are similar to [27,28]; I retain the 
modification of the gluon Fourier transform that was discussed above. ^ 

We make the following redefinition of the momentum space 1-component fields: 



The notation is as follows. In the definition of the set of 4- vectors /C, powers indicate 
how many times a or 1 appears. Underlining indicates that all permutations of 
entries are to be included. Note that the 16 lattice poles described in §2. II above are 
located in momentum space at k E (7r/a)/C. The map tt is a projection to the reduced 

•^There is no difficulty applying Reisz's power-counting theorem to pure YM. It is the SF propa- 
gator and vertices that pose problems. I include the gluon propagator (in Feynman gauge) for the 
purpose of illustrating how the YM sector is treated in tandem with SFs when PC is attempted. 
Treatment of ghosts is identical. A detailed analysis of pure YM interactions will not be required in 
what follows. 

^There is also a position space taste basis [29], which I will not discuss here, except briefly in 
For a more detailed exposition, consult also [30] as well as the reviews [31,32]. 




a 



(2.10) 
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first Brillouin zone I32a- Feynman vertices and propagators involving the fermions 
are then translated from the ICB using this identification. In practice it is helpful to 
extend the definition as follows: 

m = XAik'), m = XAik'), k = k' + -A, V k,k'eR\Ae}C. (2.11) 

a 

Here, the periodically extended definitions [cf. ()2.5|l ] of x{k),x{k) are used. 

Taking these redefinitions into account, corresponding to ()2.8p and ()2.9|) we have 
the Feynman rules [27,28]: 



= p' + ^A, q = q' + ^B)^ S'Up', q') 



= S'^W + q')— ' (2.12) 

V^^ {p = p' + ^A, q = q' + -B- k) ^ V;%(p', q'- k) 

= —giT^trW + q' + k) (e^"^ + e-'^P'^^^-^'^) ^^)a,b. (2.13) 

As in fj2.11|) . p', g', k take any values in R^. An equivalence has been used to obtain 
momentum conserving (mod 27r/a) ^Q-functions. It is reviewed in ^ and plays a 
crucial role in the generalized proof of §S1 The 16 x 16 momentum space spin-taste 
matrices (1 1)a,b and (7^ ® ^)a,b are written in the notation of [27,28]. Definitions 
can be found therein; we will not need their explicit form in what follows. 

The lattice perturbation theory also contains quark-multigluon vertices that are 
irrelevant operators, suppressed by explicit powers of the lattice spacing. They are 
important to take into account for the cancellation of divergences (see for example 
§14.2 of [31]), and are easily incorporated into the formalism that is discussed below. ^ 

3 Review of Reisz lattice power-counting 

3.1 The lattice UV degree 

A given Feynman integral is written in the general form 



,4, ,4, Vik,q;m,a) 

-. (3.1) 

C{k,q;m,a) 



thank David Adams for an important discussion regarding the role of the irrelevant vertices. 
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Here, /ci, . . . , fc^ are loop momenta and qi, . . . ,qE are external momenta. Note that 
the loop momenta are integrated over Ba^ = 'x^Ba = {—T^/a.Tt/af^; cf. fl2.2j) . Also, 
m stands collectively for mass parameters. The numerator V incorporates vertex 
factors and the numerators of propagators; C consists of a product of propagator 
denominators. 

Reisz defines the UV degree of V and C, and thence of the integral /. At higher 
orders, this must be done over Zimmermann subspaces H. To each propagator cor- 
responds a line momentum ii{k,q) (cf. ()3.12|) below). There is a sense in which L of 
these form a basis w.r.t. ki, . . . , k^, as will be explained in ^3.31 below (cf. condition 
L2). We decompose this set as follows: 

vi = ij„ VL-d = ^ji^_,- (3.2) 

We regard vi, . . . , vi-d and qi, . . . ,qE as external momenta to the Zimmermann sub- 
space H. The Ui, . . . ,Ud are regarded as internal momenta that provide a parameter- 
ization of H. We denote the set of all Zimmermann subspaces by Ti. 

The UV degree of V w.r.t. H is just the A oo scaling exponent ofVasu^Xu 
and a ^ a/X. First we define: 

V{u, V, g; m, a) = V{k{u, v, q), q; m, a), (3.3) 

using the fact that the line momenta in ()3.2|) form a basis w.r.t. k. Then, as A oo 
we extract the leading exponent:^ 

V{Xu,v,q;m,a/X) = X''A{u,v,q;m,a) + OiX"'^) <^ degr^^V = u. (3.4) 

The UV degree of C{k, q; m, a) is defined similarly. Combining the two, we have 

degTfjI = 4d + degr^V — degr^C, (3.5) 

where we recall that there are d momenta internal to H. 

3.2 Reisz's theorem 

Convergence of the Feynman integral is then proven, provided 

degF^J < 0, V Hen. (3.6) 

^Reisz adds a "hat" to the subscript of the degree operator in order to distinguish it from "the 
old polynomial degree." 
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In this case, one obtains the remarkable result: 



limi- Td^^k^^^l^ (3 7) 



where 



P{k,q;m) = \imV{k,q]m,a), E{k,q]m) = \imC{k,q;m,a) (3.8) 

a— >0 a^O 

are just the continuum limits of the numerator and denominator resp. Next I consider 
the conditions that are assumed to hold in the course of proving this result. 



3.3 The Reisz-Liischer conditions 

In the proof of ()3.5j) and its consequences, Reisz makes some assumptions about the 
Feynman integrand. The lattice power-counting theorem of Reisz has been reviewed 
by Liischer [6], and I will make use of his enumeration of the conditions that are 
assumed in the course of the proof. I refer to these as the Reisz-Liischer conditions. 

First, V satisfies: 
VI. There is an integers and function F such that 



V{k,q;m,a) = a F{ka,qa;ma). (3.9) 

Moreover, F is smooth, Ire-periodic in the momenta ka, and a polynomial in the 
masses ma. 

V2. V has a continuum limit, in the sense that 



P{k, q; m) = lim V{k, q; m, a) (3.10) 

exists. 

Liischer notes that VI — V2 are "not very restrictive." We will find that they are 
satisfied for SF-QCD. 

As stated above, the denominator function C that appears in ()3.H) is a product 
of the denominators of propagators, Ci, . . . , Cf. 

I 

C{k, g; m, a) = Y\ Ci{ii; m, a). (3-11) 

1=1 

Here, each Ci depends on a line momentum, described in more detail below. In 
SF-QCD, Ci is a trigonometric function of the line momentum £j. 
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Reisz requires that the hne momenta be natural, Defn. 3.1 in [10]. That is:^ 
LI*. The line momenta are of the form 

L E 

i,{k, q) = Y, Cr,k, + = K,{k) + Q,{q). (3.12) 

Moreover, it is assumed that Cij G Z, G R, and 

rankC,, = L, (C^, • • • , Q^) ^ V z = l,...,/. (3.13) 



L2. Define the set 

C={k^,...,kL,K,,...,Kj}, (3.14) 

where Ki = J2j Cijkj were defined in / IJ. iii)) . Suppose Ui, . . . ,ul are linearly indepen- 
dent elements contained in C. Then the loop momenta can be expressed in terms of 
them with integer coefficients: 

L 

ki = 'YAijUj, AjEZ. (3.15) 

i=i 

Note that this property was used above in ()3.3p . It is in this sense that the hne 
momenta appearing in ()3.2j) form a basis w.r.t. ki, . . . , ki- 

The foUowing requirements are imposed on the functions that appear in ()3.11|) : 

CI. The propagator denominators can be expressed as 

Ci{£i; m, a) = a^'^Gi^aif, am), (3.16) 

where Gi is a smooth function that is 2-^ -periodic in the momentum iia. Also, Gi is 
a polynomial in the mass ma. 

C2. The denominators have the conventional continuum limit: 

hm Ci{^u a) = + m\, (3.17) 

for fixed li,mi. Here, mf is a homogeneous quadratic polynomial in the mass param- 
eters of the theory. 

^The asterisk denotes that I have modified Liischer's statement of the condition in order to bring 
it into harmony with stricter definition given by Reisz. 
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C3. There exists an > and an A > such that 

IQiiu m,a)\>A + m^) , = J2^ sm\ii^a/2), (3.18) 

for all a < and all ii E Ba- 
in Reisz's proof, he does not require C3, but instead a condition that has the 
same effect — a lower bound on the propagator denominator: 
C3*. There exists an aQ > and an A > such that 

\Q{if,m,a)\>A{e^^+m^,), (3.19) 

for all a < ao and all £i E Ba- 

Note that the r.h.s. of the inequahty is a continuum expression. The reason that 
this is equivalent to C3 is that C = if + mf itself satisfies C3*. Thus we can always 
replace the bound in C3 by the continuum (rational) expression in C3*. In fact, the 
essence of Reisz's proof is to replace lattice expressions by bounds that are rational 
and have a continuum interpretation. For this reason I prefer C3*. 



4 The conditions of Reisz's theorem vs. staggered 
fermions 

Here I examine the Reisz-Liischer conditions in relation to SFs. It will turn out that 
the Reisz-Liischer conditions fail in both the ICB and the MSTB. The essential reason 
is a mismatch between the Kahler-Dirac lattice 2aZ'^ and the gauge lattice aTl^ . (The 
former is the natural lattice on which to formulate free SFs, whereas the latter is the 
lattice on which the pure YM theory is formulated.) 

In these considerations, it is implied that the 8 functions that appear in the Fey- 
man rules of §^2312121 have been integrated against (except for the overall 5 function 
that always occurs), leading to the line momenta ^i(fc, q). 



4.1 1-component basis 



Consider the denominators of the propagators D'^^{k,i) and S^^{p,q) in relation to 



the Reisz-Liischer conditions. For the gluon (B) and quark (F), 



C 



B 



-2 



5^4sin^^ + (Ac 



-2 



a) + (ma) 



(4.1) 
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where ^ is the hne momentum flowing into the propagator. Both are of the form 
and are periodic on Ba] i.e., unchanged under £ ^ £ + (27r/a). Thus CI is 
satisfied. It is also obvious that C2 holds. Whereas C3 holds for Cb-, it does not 
hold for Cp. This is because the latter has lattice poles away from the origin of Ba, 
as I now show. 

The proof consists of showing that there exist ^ & Ba such that C3 fails. In 
particular, suppose that £ = {TT/a)B, where B E K,* = K, — (0^). Then Cp = 'm? and 
= 4a~^^^_B^ = 4a~^|i?|, where l"^ was defined in ()3.18j) . Let A be any strictly 
positive real number. To satisfy C3 it is necessary that for sufficiently small Oq, and 

any a < aQ 

- A{4:a-^\B\+m^) >0. (4.2) 
This can only be true for A < 1. But, for any a such that 

AA\B\ 



< a < a' 



_{1-A)m^ 



1/2 

(4.3) 



condition ()4.2j) is violated. Thus we can never choose ao small enough to satisfy C3. 
Put simply, near one of the extra lattice poles, Cp = 0{m'^) whereas A{i'^ + m^) = 
74(9(a~^); so, for small enough a the latter is always larger. 

The numerator of the quark propagator ()2.8j) will contribute to the Feynman 
numerator V in (jH.lj) . and it is easy to see that it satifies V1-V2. The (5a-functions 
that appear are periodic on Ba by construction. It is always possible to choose 
loop momenta ki and line momenta ii such that p = ii{k,q) in the denominator 
of ()2.8j) . where q is external momentum passing through the propagator. Due to 
the vr/a violations of momentum in the ICB, it is not guaranteed that the ii are 
natural; it is therefore possible that L1-L2 are also not satisfied. In any case ii = 
Cijkj + Qi{q) + Zi{n/a) with Cij G Z and Zi G Z^. Then under kj — > kj + {2t[ /a)v, for 
any loop momentum kj, the numerator term ^^sin£j^a is unchanged. The vertex 
V^'^^{k;p,q) also satisfies VI, because it is has been constructed to be periodic on 
Ba- It is obvious that the vertex satisfies V2. 

Thus we see that the Reisz-Liischer conditions fail to hold in the ICB principally 
for the reason that has been pointed out by Liischer [6]: the fermion propagator has 
too many poles in Ba- Also worrisome is the vr/a violations of momentum. The latter 
problem will be eliminated in the basis that I discuss next. 
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4.2 Momentum space taste basis 

Note that the denominators of the MSTB propagators are the same as in the ICB. 
One has exactly the same violation of the Reisz-Liischer conditions as in the ICB, 
due to additional quark poles in Ba- 

One might think to instead apply the Reisz-Liischer conditions on the reduced 
Brillouin zone B2a, defined in (I2.1(J|) . This effectively replaces the lattice spacing a 
hj b = 2a. Whereas C3 is satisfied if Ba is replaced by Bb = B2a, the Feynman rules 
do not enjoy 27r/6 = vr/a periodicity. Thus VI and CI would be violated if we took 
this approach. In fact, in the generalization that is introduced in §H1 below, the lack 
of 7r/a-periodicity will be addressed "head-on". 

4.3 The mismatch 

The problem with 7r/a-periodicity is an inevitable consequence of the mismatch be- 
tween the Kahler-Dirac lattice and the gauge lattice. To see this, note that the 
Fourier transform ()2.4|) has been formulated w.r.t. the translation invariance group 
of the gauge lattice, generated by shifts fia. As a consequence, the fields (^(fc) are 
periodic on the reciprocal lattice {27!-/a)Z'^. If not for the gauge fields, we could per- 
form a Fourier transform w.r.t. the translation invariance group of the Kahler-Dirac 
lattice, generated by shifts fib = 2jla. To accomplish this, we pass to the position 
space hypercube basis [28-30] before taking the Fourier transform: 

X(r = 2y + A)= XA{y). yeZ\ AeK,- 

XAiy) = [ ^.e^'^'xAik) = [ ^.e^'^^'-'xAik). (4.4) 

It is easy to see that the quark propagator for the Xa(^) fields is vr/a-periodic. This 
just follows from the inverse Fourier transform: 

XA{k) = b'Yl ^""'"xAiy), (4.5) 

which clearly leads to 

XA{k + i7r/a)z) = XA{k) V z e Z\ (4.6) 

Thus the Feynman rules for XA{k + {j: la)z) and XAik) will be identical. This is just 
to say that the reciprocal lattice of the Kahler-Dirac lattice is (27r/6)Z^ = {n/a)7J^. 
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5 The Reisz proof generalized 



I now extend the Reisz power-counting theorem and proof such that the deviations 
from the Reisz-Liischer conditions can be overcome. The trick is to use: 

(i) the MSTB rules with momentum-conserving 5-functions, 

(ii) integration extended to full periods, 

(iii) resolution of identity on the reduced reciprocal lattice, and 

(iv) tranformed Feynman rules that absorb loop momentum shifts 



The Feynman integral is written in the MSTB. The domain of integration is Ba, 
using the trick of The vertices and propagators conserve momentum, using the 
equivalence of ^ Of course, one first integrates and sums against all 6 functions. 
This should be done in such a way as to conserve momentum as it flows through the 
diagram, using the 27r/a-periodicity wherever necessary.*^ Then we are left with just 
integrals over loop momenta, which can be routed such that the line momenta are 
natural. That is, the momenta are routed just as in a continuum Feynman integral. 

5.1 The denominator 

1 now describe the modiflcation that overcomes the principal difficulties posed by SFs. 
One breaks up the line momenta into those corresponding to bosons (gluons) and 
fermions (quarks): if , . . . , if^ and if , . . . , if^ resp. Then one inserts into the Feyn- 
man integral, for each ii, the resolutions of identity that are described in Isiif) 
deflned in (IA.2|) and Ipiif) deflned in ()A.4|) . One obtains an expression analogous 
to Reisz's Eq. (4.4) [10] — a sum of integrals that comprises a domain decomposition: 



^This is optional because the MSTB Feynman rules only conserve momentum mod l-nja. 



ki ki + {7T/a)Ai, Ai G /C. 



(5.1) 




Jb<^{1,...,Ib}, JfC{1,...,/p}, 

zb = {zBi\i e Jb), zf = {zpjlj e Jf), 



(5.2) 
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with individual terms of the form: 

d'^^k V{k,q]m,a) 



B.- nfil CB(^f ; A, a) Cp{q- m, a) 



xn0(r-ii^f-^^-^n)n0f(^f)- (5-3) 

Note that J collectively denotes J^? and so on. The decomposition has the fol- 
lowing intuitive meaning: £j G J are "e-near" to a lattice pole, whereas li ^ J are 
"e-far" from a lattice pole. 

For e, a sufficiently small, the arguments of Reisz's Appendix D [10] extend in an 
obvious way to show that there exists /c^^-* = {kf \ . . . , k^^^) G BJ" s.t.: 

Kf{k(^)) = —ZB^, Kfik^'y) = -ZF„ ieJb, jeJp. (5.4) 
a a 

Note that Ki{k) was defined in ()3.12p . Using the fact that are natural, it is a trivial 
extension of Reisz's Lemma D.2 [10] to prove that there exist reduced reciprocal lattice 
vectors 

Ai,...,A^G-Z^ (5.5) 

ct 

such that for i E Jb, j & Jp 

KfiA) = —ZB^, Kf{A) = -ZFy (5.6) 
a a 

The Aj are determined in terms of a basis chosen from {Kf ,Kj}, as explained in 
Reisz's Appendix D [10]. Thus we define new loop momenta k[ through: 

k, = k[ + /\i = k'i + -5i V i = l,...,L, (5.7) 
a 

where in the last step integer-valued 4-vectors 5i have been introduced for future 
convenience, following Reisz. A new domain of integration results: 

dj = |A;' G R^^ - - - A,^ < k[^ <-- aA , (5.8) 

identical to Reisz's Eq. (4.7) [10]. 

For the line momenta £f G JbA^ ^ Jf, ()5.7j) has the effect 

(fc) = (^') + -zf. if{k) = if{k') + -zf, (5.9) 
a a 
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a generalization of Reisz's (4.5) [10]. When this is accounted for in ()5.3p . the Heaviside 
step functions in ()5.H|1 just force if{k') G Jb, ^f{k') G Jp into the e-neighborhood of 
the (unique) pole in Ba and respectively. As a consequence the following bounds 
hold: 

Cp\if e Jf) < apiifik'f + mY\ (5.10) 

generalizations of Reisz's (4.8) [10]. Here, aBjOp are constants that always exist for 
e, a sufficiently small. For ii J, the line momenta are outside of the balls of radius 
en /a that are centered on sites of the (reduced) reciprocal lattice for (quarks) gluons. 
Therefore they are bounded by: 

C-i(£f ^ Jb) < iBa\ C^\£f ^ Jp) < lFa\ (5.11) 

generalizations of Reisz's (4.9) [10]. Here, ■Jb,1f are constants that always exist for 
e, a sufficiently small. 

For the line momenta if ^ JB,if ^ Jp, the shift ()5.7|1 is only guaranteed to have 

£fik) - if{k') = C^A„ G -Z^ if{k) - if{k') = Cf„A„ G -Z\ (5.12) 

a ■'a 

Whereas 6f is vr/a-periodic, the function 6f is only 27r/a-periodic. Some explicit 
dependence on A will result, and will be addressed below. Gathering together the 
various results, we can bound ()5.H|1 by: 

ijz < a1fa''/{^Ba''f^-''^\iFa'f^-''^^ 

X 



f d^^k' V(A;' + A,g;m,a) 

J (J J 

n [{(if{k'f+xr'Q G^-ii^f(^') 

n kW + m2)-i0(^e-||£f(A:')|| 



X 



X 

X n + Cf„A„) n 0f (^f (^'))- (5.13) 

Here, the hBihp are the number of elements in Jb,Jf resp.; i.e., the number of 
line momenta that are e-near to lattice poles. The denominator has been expressed 
entirely in terms of rational functions. The numerator V requires further study: the 
shifted loop momentum argument k' + A can be accomodated into Reisz's techniques 
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to bound the numerator, as will be discussed further in ^5.21 and ^5.41 below. It will 
be seen that the 9, 9f , 9f functions do not pose any difficulty, as they just restrict 
the domain of loop integration. With the bound in the form ()5.13|) . it is quite simple 
to extend the remainder of Reisz's manipulations. Using them, I will formulate and 
prove the SF power-counting theorem. 

5.2 The numerator 

The magnitude of the SF numerator is also easy to estimate, using the decomposition 
Jp and the shifted line momenta ij{k'). However, cancellations associated with the 
spin-taste algebra will be important to take into account in order to get the correct UV 
degree for a given diagram. For this reason it is better to abide by Reisz's approach 
and treat the numerator V whole. 

I now make a few remarks regarding the effect of the shift A that appears in the 
numerator of (j5.13|) . This will lead to modifications of propagator numerators and 
of vertex factors. In the words of Reisz, a generic shift Aj G "would produce 
explicit negative powers in the lattice spacing destroying convergence." However, 
though the shift involved here is not an invariance of the Feynman integrand, it is 
nevertheless special: Aj G (7r/a)Z^. As I now discuss, it is possible to eliminate this 
explicit through a transformation in the form of the Feynman integrand. The 
transformed integrand trades sin ^ cos in various places, and/or introduces factors 
of (—1). Furthermore, the number of possibilities for how the propagators and vertex 
factors are transformed is finite. The UV degree of the transformed numerator is 
then determined in accordance with Reisz's definition. This degree is then used in a 
generalized computation of the UV degree of the integral J, as will be seen in ^5.31 
below. 

As an example of the transformation induced by A, consider the quark-gluon 
vertex ()2.13|) . For the sake of argument, suppose that each line entering the vertex is 
internal, with 

p^ii, q ^ i2, k ^ is, 

ii = Ci,k, + Qi{q) V « = 1,2,3. (5.14) 
Thus in the redefinition ()5.7|) the momentum-dependent factor in the vertex trans- 
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forms as: 



(5.15) 



Note that since the hne momenta are natural, Cij G Z. It follows that, as promised, 
factors of (—1) have been introduced. Momentum conservation implies 

Cu + C2i + = 0, (5.16) 

but this still allows for the transformation ()5.15|) to have a notrivial effect on the 
quark-gluon vertex. In the case that C^i5i^ = 1 mod 2, the factor cos{ii^{k)a + 
(l/2)£3^(A;)a) is exchanged for sm{ii^{k')a + {ll2)l^^{k')a). [Here, it is implicit that 
overall exponentials are factored out to rewrite the expression using trigonometric 
functions.] Since the latter starts at 0{a), rather than 0(1), the UV properties of 
the vertex are changed in a significant way (lowered by 1). This does not destroy the 
convergence of the numerator; in fact, it improves it. 

As another example, consider the triple-gluon vertex. In the conventions taken 
here,^ 

Kp^(P' 5' ^) = f /"^'V^a(p + q + k) {e^^^ - e-'<^^+''^^) (1 + e-^'^(P'^+^^)) 

+ cyclic permutations. (5-17) 

Applying ()5.14j) . one finds: 



iat-j,p(k) _ -ia{hp{k)+i2p{k))\ /i I „~ia{h^,{k)+£3t,{k)) 



)(l + 

_^C3i5ip {^ial3p{k') _ ^_-^C2iSip^-ia{e3p(k')+e2p(k')) 



e 



(5.18) 

Here again, the explicit contained in Aj is traded for factors of (—1), that in some 

cases interchange sin ^ cos. Raising or lowering the UV degree by 1. 

An increase in the UV degree of the numerator under transformations such as 

(|5.18p will not "destroy" the "convergence" of the Feynman integral. Rather, it will 

only make manifest the cutoff dependence. To the extent that a Feynman integral has 

positive UV degree, subtractions are required in order to have a convergent result, 

regardless of the basis of loop momenta. 

^See for example Eq. (3.225) of [32], taking into account the 27r/a-periodic Fourier transform 
convention H2.4|l that I use. 
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Quite generally, the factors in the numerator V{k, q; m, a) are trigonometric func- 
tions that are 27r/a-periodic. The redefinition ki = k'^ + {7c/a)6i, 6i G always results 
in a half- or full-period shift. The rule is that, prior to computing the UV degree, one 
should eliminate the explicit vr/a factor using elementary trigonometric identities, as 
has just been illustrated for the quark-gluon and triple-gluon vertices. 

The best strategy to deal with this is to extend the Feynman rules to incorporate 
n/a shifts. Then the integral (j5.13|) should be interpreted in terms of these new rules. 
Symbolically, 

V{k' + (7r/a)(5, q; m, a) = V{k' + (7r/a)y4, g; m, a) = VA{k', q; m, a), 

= Si^ mod 2, Ai e /C. (5.19) 

Here, Va is written in terms of the generalized Feynman rules. (Note that in the 
second step the 27r/a-periodicity has been used to express the numerator V in terms 
of the transformed one with an index restricted to A G IC^.) Once this has been done, 
all of Reisz's techniques for the UV degree analysis of the numerator apply. This can 
be seen from the fact that the numerator satisfies the Reisz-Liischer conditions, after 
the explicit factors of vr/a have been eliminated. The essential reason for this is that 
sine and cosine are analytic functions. 

5.3 The generalized theorem 

These considerations lead to the following generalization of Reisz's theorem: 
Definition. Let = Va/Ca, A G /C^ denote the transformed Feynman integrand. 

That is: 

VA^k, g; m, a) = V{k + {n/ajA, q; m, a), 

CA^k, q; m, a) = C{k + {j[ /a)A, q; m, a). (5.20) 
Generalize the UV degree as follows: 

degr^F = max degr^F^, degr^F^ = degr^VA - degr^C^; 

degr^J = 4d + degr^F. (5-21) 

Recall that ui, . . . ,Ud parameterizes the Zimmermann subspace H. 
Proposition. Suppose that 
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degr^J<0 V H eU. (5.22) 

Then / converges, and 



oo 



,„„/= I d"-k^^^^^, (5.23) 



where 



PA{k,q;m) = \imVA{k,q;m,a), g; m) = hm C^(A;, g; m, a) (5.24) 

a^O a— »0 

are just the continuum hmits of the numerator and denominator resp. This indicates 
that various regions of loop momenta may contribute to the continuum hmit, due to 
the presence of doublers in the fermion spectrum. 

5.4 Proof 

Starting with ()5.3p . one makes the redefinition ()5.7p . Then the numerator is replaced 
by VA{k', q; m, a), as in ()5.19p . Once this has been done, Ijz is in the form considered 
by Reisz. Due to the assumption ()5.22|1 . the remainder Ra in the decomposition 

Va(A;, q; m, a) = PA{k, q; m) + RA^k, q; m, a) (5.25) 

does not contribute in the continuum limit, as follows from Reisz's arguments in §7 
of [10]. Thus one can replace Va by the rational function Pa in the numerator of 
ijz- Furthermore, Reisz's arguments show that the Ijz term that maps to the index 
Ae JC^ just yields 

rf-fc ^4^lii4 (5.26) 
EA{k,q;m) 

in the continuum limit. The result ()5.23p follows immediately. 

To clarify this, I discuss some of the details of §7 of Reisz [10]. Taking into 
account the decomposition (j5.25p and the effect of (j5.7p and (j5.19|) on (j5.3|) . we have 
the decomposition 

/j. = 4 + jji, (5.27) 



where 



Jz 



d^^k' PA{k',q;m) 
aj CA{k', q; m, a) 

xlle{le-\\ifik')\\) n0(^^-ll^f(^')ll 

X n ^fi^fif^') + CL^rn) n 0f «(^')), (5.28) 
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and Jj^ is defined with P4 replaced by Ra- Note that there is a correspondence 
between the index A and the index Jz, since the latter determines the shift ()5.7j) that 
is required. When one compares to Reisz's (7.1)-(7.2) of [10], the only difference that 
one finds is in the factor rii^js ©f (^f (^') + Q^^m)) with its extra argument Cj^A^- 
This only restricts the domain of integration; it can only make the integral smaller, 
so Reisz's bounds still hold. E.g., 

^ f d^'k' \PAik',q;m)\ 

where kj is the same domain as in Reisz's (7.6) [10]. (Following Reisz, the domain 
()5.8|) is extended, cxj ^ k,j, which can only increase the value of the bounding contin- 
uum integral Ij^-) Similar remarks apply to If^. Note that this is just an exploitation 
of the bound (jHlSl)- 

It follows that the remainder of Reisz's arguments of §7 [10] apply, with only the 
following modification. It is still true that Ij^ has a nonvanishing continuum limit 
only if all line momenta are near poles: Jb = {1, . . . , Ib} and Jp = {1, . . . , If}- 
However, rather than Reisz's Eq. (7.18) [10], we must write {I = Ib + If and i 
collectively denotes all line momenta): 

L 

Ki{k) = V Ci,k^ = -Zi, t = l,...,I, (5.30) 
^-^ a 

i=i 

which differs in that 27r/a has been replaced by vr/a for the coefficient of Zi. Just as 
in his case, this equation has a unique solution, due to rank(Cij) = m. Due to the 
naturalness of the line momenta, the solution is of the form 

h = (vr/a)5,, 6, e Z\ (5.31) 

These are nothing but the 6i that are used in the shift (j5.7j) . 

Now consider the possible solutions to (jSISH), given ki G Ba- Unlike what occurs 
in the case considered by Reisz, there are multiple possibilities for 6i that will work. 
In fact, they are nothing other than the G /C. That is, the sum over AelC^ that 
is taken in (j5.23p is in one-to-one correspondence with the set of fcj = {71/0)81 that are 
solutions to (j5.31|) . (As will be seen in the example of ^5.5.21 it can happen that for 
certain A = 5 G /C^, we get Zi 2Z^ where i corresponds to one of the gluon lines, 
C.f . In that case the contribution I a will vanish, since the gluon is always far from 
the pole in /j^.) From this we see that the sum in ()5.23|) is just the modification that 
is required to extend the arguments of Reisz's proof. 
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5.5 Simple examples 

5.5.1 1-loop with external gluons 

Consider any 1-fermion-Ioop diagram with external gluon legs.^° For example, the 
4-gluon diagram Fig. [T] Recall that in the MSTB rules, a factor of 1/16 is supplied 
for each fermion loop, due to the extension of the integration domain. On the other 
hand, when one takes the continuum limit of the numerator and denominator, only 
one pole region {k ^ 0) occurs in the continuum Feynman integrand [Eq. ()5.26p with 
A = (0,0,0,0).]. It is only when we include all 16 contributions that come from the 
sum over A & JC that appears in ()5.2H|1 that we get the correct overall factor. 

In detail, the denominator of the integrand is invariant under the shift of loop 
momentum, k —>■ k + {7i/a)A, since it involves only the SF propagator. Taking into 
account the change in the SF propagator numerators and the gluon-quark vertices, it 
can be shown that the change in the numerator of the integrand is equivalent to just: 

(^¥T)--(-)^''(^^), (5.32) 

in every place that (7^ ® 1) appears, and everything else left as it was before the shift 
in loop momentum. This is just an automorphism of the Clifford algebra, and will 
not change the value of the character — i.e., the trace that occurs in the numerator. 
Thus each of the 16 Feynman integrands, A G /C, are identical in their continuum 
limit to the one obtained at A = 0. This just cancels the factor of 1/16 to give the 
desired result. 

5.5.2 1-loop with external fermions 

Next consider diagrams such as Fig. |21 In this case the gluons contribute to the 
denominator. When the shift k k + [it /a) A is performed, with A 0, the gluon 
denominators are moved far from their poles when /c ~ 0. It follows that C^^ ~ a^, 
where Ca is the transformed denominator of the Feynman integrand. The numerator 
of the integrand Va is well-behaved as a ^ 0. Thus for A ^ the integrand vanishes 
in the a — limit. This is just as required: there is no 1/16 factor to be compensated 
because we do not have a fermion loop. Only the A = contribution survives in the 
a —>■ limit, and the result is just the usual continuum expression. 

Vacuum polarization will of course require a subtraction for the power-counting theorem to hold. 
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Figure 1: Example of finite diagram where all 16 choices of ^4 e /C contribute. 



Figure 2: Example of finite diagram where only 1 choice of ^4 e /C contributes. 
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6 Discussion 



In this article I have demonstrated how to extend the techniques of Reisz's power- 
counting theorem to the case of staggered fermions. It is fortunate that with a few 
straightforward modifications, the bulk of Reisz's arguments apply. It is of great prac- 
tical importance that his power-counting theorem for generalized continuum Feynman 
integrals (§5 of [10]) continues to be applicable. One thing that remains to be done 
is to use the staggered fermion power-counting theorem to prove perturbative renor- 
malizability, following [7,8]. Also, applications of the theorem to higher orders in 
perturbation theory should be explored in more detail. Generalizations to other sorts 
of theories that contain doublers could also be considered, since the technique that 
has been introduced here is not very specialized. It is worth emphasizing that all of 
the manipulations that were performed here apply equally well to improved staggered 
fermion QCD. 
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Appendices 



A Resolutions of identity 



A.l Resolution on Ba 

In §4 of [10], Rcisz introduces the following step function: 



if [[£ - ^^11 < for some z e 

1 otherwise. 



(A.l) 



A superscript B has been affixed to distinguish it from another step function that 
will be defined below. For any i one can resolve identity as: 

2n 



i = iB{i) ^ ef(£) + 5]ef-6- i--z) 



(A.2) 



Here, O is Heaviside's unit step function. 
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A. 2 Resolution on 

Define a step-function analogous to Reisz's: 



ef(£) 



if ||£- ^zll < for some z G Z^, 

1 otlierwise. 



For any i one can resolve identity as: 

1 = 1^^) = 0fw + EQ(^^ 



TT 

z 

a 



(A.3) 



(A.4) 



As above, G is Heaviside's unit step function. This resolution is useful for line mo- 
menta of SFs, since it isolates the regions that are near SF poles. 



B Simplification of MSTB rules 

Here I establish a very important simplification that is used to derive the MSTB rules 
()2.12|) - ()2.13|) . It was employed, for instance, by Patel and Sharpe [28]. 



B.l ^-function transformation 

Without loss of generality, for any momenta p,q & we can write 

p = p' + -A, q = q' + -B, p',q'eR\ A,BelC. (B.l) 
a a 

That is, p', q' are not necessarily in the reduced first Brillouin zone B2a- Even if p, q 
are integrated over we can impose that p', q' are also integrated over Ba by using 
the equivalence that is established in ^jCl In fact, this is exactly what I do when 
passing to the MSTB from the ICB. It is also what has been done by Patel and 
Sharpe [28]. 

It is not hard to check that the 5-functions that appear in the ICB quark propa- 
gator ()2.8p and vertex ()2.9|) obey the identities 

5a{v' + q' + k + ^A + B + Jj)) = 

4 

n P-^IP'. + ^^ + k.)5l+B.+,.fl + 4 {v'. + + + 3 , (B.2) 

v=l 
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with A; = in the quark propagator ()2.8|) . and (5'^' a Kronecker S mod 2. It is worth 
noting that in the transition from ^^-^^ to ^U^ - ^H^ . the identity^^ 

4 

is used. The factor (— )^'' is obtained from the p-dependent prefactors of the 6a- 
functions in p.8|l - p.9|) under p = p' + (7r/a)A. However, the other (5-functions that 
occur in ()B.2|1 need to be taken into account. Next I will demonstrate the equivalence 
that allows us to eliminate the 6 functions that violate momentum conservation by 
vr/a. In the proof the power-counting theorem, this is key to obtaining natural line 
momenta in the continuum, bounding integrals Ij^ and if^. 

B.2 Equivalence 

As mentioned above, the trick of ^is used to extend the integration of p', q' to Ba- 
Then a typical integral against the 5a{- ■ ■ + vr/a) parts of (jB.2j) takes the form: 




where / represents the rest of the integrand, where ". . ." corresponds to integrals over 
other momenta components and other momenta. I hide integrations over these other 
variables, for simplicity of notation. The functional form of /, namely the way that 
p' and A appear together, etc., is guaranteed by the fact that we start from the ICB. 
Next I make the redefinitions 

ql-^ql--, B,^B, + 1, (B.5) 
ct 

to obtain 




x4 (pI + ql + K) 521+B.+/i.,o / {p'u + ^A,, q'^ + . . .) . (B.6) 

^^This identity is easily checked using, for example, the definitions given in §2.1-2.2 of [28]. 
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Finally, since the integration over is a full period, we can shift the domain 

/ - / dq'^ (B.7) 

JQ J-Tx/a 

at no cost. This establishes the equivalence 5a{- ■ ■ + vr/a) to the momentum conserving 
terms 5a{- ■ ■)■ In essence, we trade momentum violation by 7r/a for a shift in taste. 
Thus we find that under the integration d'^q' , 

4 

n [^-^P'u + <i'u + K)5t+BM + (P'u + + + 3 ^AM+f..,l_ 

i&Up' + + ^)4'U+M,o- (B-8) 

We could also have used the p' integration in these manipulations. If 2 momenta 
are external, we cannot use this equivalence and more care is required [33]. This 
technicality does not affect the power-counting considerations here, because we are 
only interested in IPI loop diagrams (the aim is to study renormalization) , which do 
not contain vertices with 2 external momenta. 



C Domain extension in MSTB 

Here I prove that after the decomposition (jB.lll of momenta for the transition ICB 
— >■ MSTB, we can double the domain of integration for p', etc. The point is to obtain 
integration over a full period. It also is invoked in the simplification of 

In the manipulations of this paragraph, p will denote any of the components 
of a momentum that is integrated over in the Feynman rules in the ICB. We use 
equivalence of integration over any full period to write: 




(C.l) 



where in the last step the integration has been prepared for further manipulations. 
Then we decompose pi,P2 onto separate reduced domains: 

Pi = p[ + -Ai, p[ e (-7r/2a, 7T/2a], A, G {0, 1}, 
a 

p^=p'^--A2, G (7r/2a,37r/2a], ^2 G {0, 1}. (C.2) 
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The integration of ()C1|) can then be expressed equivalently as: 

1 / /"T/Sa /'37r/2a \ -r /■37r/2a -r /"Tr/a 

^ \J-K/2a J-K/2a / ^ ^ J-7r/2a ^ ^ J -w/a 

where in the last step I have again used the equivalence of integrations over a full 
period. 

Extending this manipulation to all components, we obtain the identity 

/ ^V---) = ^E/ ^V(---W'+fA. (C.4) 
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